Parity effect in a small superconducting grain: A rigorous result 
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The parity effect in an ultra-small superconducting grain is examined. By applying a generalized 
version of Lieb's spin-reflection positivity technique, we show rigorously that the parity parameter 
Ap is nonvanishing in such a system. A positive lower bound for Ap is derived. 
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I. INTRODUCTION 

Recent experiments by Ralph, Black and Tinkham 
jl],^| on tunneling currents through nanometer-scale Al 
grains have rekindled theoretical interest in modifications 
of the original Bardeen-Cooper-Schrieffer (BCS) theory 
H when the average level spacing Se becomes comparable 
to the bulk superconducting gap A. Although the pos- 
sible failure of the BCS mean-field theory in this regime 
was speculated by Anderson in 1959, quantitative 
studies of the crossover from ultra-small grain to bulk 
behavior emerged only recently J^-|l^] . 

A natural generalization of the superconducting gap 
A for a bulk superconductor to the case of a small grain 
with a fixed number of electrons is the parity parameter 

A P ee E (2N + 1) - ^[E (2N) + E {2N + 2)] (1) 

where Eq(M) is the ground state energy of a grain of M 
electrons. Janko, Smith, and Ambegaokar [|| first dis- 
cussed the parity dependence of the superconducting gap 
in a small superconductor and showed that Ap reduces 
to the bulk superconducting gap A when the level spac- 
ing Se — * 0. Using a path integral technique, Matveev 
and Larkin j|] calculated Ap for the BCS Hamiltonian 
[see Eq. and found a nonvanishing result also in 

the ultra-small grain limit Se 3> A. This is in contrast 
with previous mean-field calculations |6|,Q] which yielded 
a vanishing superconducting gap when the grain becomes 
sufficiently small. Recent numerical investigations [ pd|Jl2| 
on specific models in the region Se « A have confirmed 
the conclusions of Matveev and Larkin, although a full 
analytical solution of the problem even for the simplest 
type of models in this class is not yet available. 

The purpose of this paper is to present some rigorous 
results on the positivity of the parity parameter A p un- 
der somewhat more general conditions than those con- 
sidered so far in the literature. By applying the re- 
cently developed Lieb's spin-reflection positivity tech- 
nique |l4|-|l6[] , we show that Ap is strictly positive in the 
BCS model of superconductivity for any value of Se/A 
and an arbitrary distribution of the single-particle en- 
ergy levels £fc (see below). In addition, we derive a lower 
bound for Ap which is independent of the number of 
particles in the grain. 



The paper is organized as follows. In Sec. II we intro- 
duce the model and collect some of the known symmetry 
properties of the Hamiltonian considered. In Sec. Ill we 
introduce a generalized version of Lieb's spin-reflection 
positivity method and present a proof for the strict pos- 
itivity of the parity parameter Ap. Specializing on the 
commonly studied case of a constant pair coupling con- 
stant, we derive a positive lower bound for Ap. Sec. IV 
contains a summary of our main results and discussions 
on the further possible extensions. 



II. THE HAMILTONIAN AND ITS BASIC 
SYMMETRIES 

As a plausible description of electronic states in an iso- 
lated grain of a superconductor yV(, let us consider the 
BCS pairing Hamiltonian Q 

k , " Vk\, \e k r\<ui D 

(2) 

Here k labels the single-particle eigenstates | k) of M. 
with energies and c ka (cka) is the fermion cre- 
ation (annihilation) operator which creates (annihilates) 
a fermion of spin a in state | k). The cut-off energy for 
the pairing interaction is given by a; pi > 0, which, for con- 
ventional superconductors, can be taken to be the Debye 
frequency. The sums in (§) extend over all admissible 
states | k) of M which satisfy the condition | e^ \ < lod . 
The total number of such states is denoted by Ne- The 
coupling constants g(k, k') are assumed to be positive 
for all pairs of admissible states. In the Hamiltonian (j^), 
only those matrix elements of the interaction responsible 
for superconductivity are included. The contributions of 
the other terms are assumed to be negligible. In partic- 
ular, no pair breaking terms are present. 

The average energy level spacing Se = l/A/"(ep) in a 
metallic grain of diameter d is proportional to where 
A/"(ep) is the density of states at the Fermi level. In the 
bulk limit d — > oo, the BCS mean-field theory applies 
and a superconducting ground state, characterized by a 
finite superconducting gap 



1 



A ~ 2lod exp(—Se/g), (3) 

is obtained. However, when the grain size is reduced, Se 
increases and eventually becomes comparable or larger 
than the bulk superconducting gap A, at which point 
quantum fluctuations in the occupation numbers need to 
be taken into account. 

The analysis presented in this paper is for a grain with 
a fixed number of particles. It is therefore important to 
note the following symmetries of the Hamiltonian (|^), 

(i) Hm commutes with the total fermion number op- 
erator N = J2k (™&T + wnere n ka = c\ a c ka is the 
fermion number operator on state | k) with spin a. Con- 
sequently, the total Hilbert space of the Hamiltonian de- 
composes into numerous invariant subspaces {V(M)}. 
Each of them is characterized by the particle number 
M. 

(ii) Next, consider the total spin operators of the grain 
M, 

k k 

Sz = \^{nk] -fi H ). (4) 

k 

Then, Hm also commutes with these operators. There- 
fore, both the total spin S and the z-component S z are 
good quantum numbers. In addition, for a given M and 
S, states with different S z are degenerate. 

(iii) Another important property of the Hamiltonian 
is that it does not contain pair-breaking or single- 
particle hopping terms. Thus, in constructing the eigen- 
states of the Hamiltonian, one can further restrict oneself 
to a subspace with a fixed proportion of paired and un- 
paired particles. In this subspace, the unpaired particles 
go into a fixed set of single-particle eigenstates, while the 
paired particles may choose to be in any of the remain- 
ing states. This property greatly simplifies the numerical 
procedure Jll],[l2j to find the ground state of the Hamilto- 
nian (^|), although in the analysis presented in this paper 
we do not need to presume such behavior. 

In Refs. M , [jllj and Jl2| , the case of uniform level spac- 
ing e k = kS was considered. Furthermore, g(k, k') = g > 
is assumed to be a constant for any pair of admissible 
states | k) and | k 1 ). In the following, these conditions 
will be relaxed. We shall only require that {e k } are real 
quantities and g(k, k') > for any pair of admissible 
states | k) and \ k'). 

III. THE PARITY PARAMETER 

A. Matrix representation of wavefunction 

Our result on the positivity of the parity parameter 
Ap is derived from a special representation of eigenstate 
wavefunctions previously used in the study of the bind- 
ing energy of fermions in the negative-?/ Hubbard model 



]17[ with an odd number of particles. The idea originates 
from Lieb's spin-reflection positivity method ||i^-|l6f for 
a system with equal number of up and down spins. In the 
following we repeat some of the essential steps leading to 
the representation and refer the reader to the literature 
for additional details. 

By regrouping the fermion operators, we can rewrite 
the Hamiltonian (j^) as 

Hm = ^ tkC ka Cka 

fc, <T 

9Qt, k ') (cIt^'t) ( c li c k'l) (5) 

\e k \, \e k ,\<uj D 

Notice that the Hamiltonian is symmetric with respect to 
the up-spin and the down-spin fermion operators (spin- 
reflection). However, this form is not particularly conve- 
nient when we calculate matrix elements of Hm in an oc- 
cupation number representation. To deal with this prob- 
lem, we shall consider the so-called "pseudo-fermion" 
operators, which were first introduced by Janko et al. 
|5| to study systematically BCS superconductivity with 
fixed number parity and then, re-introduced by Lieb and 
Nachtergaele jl6| to show that the ground state of the 
Hubbard model at half-filling is the global ground state 
when the chemical potential \i = ^. These pseudo- 
fermion operators are defined by 

(7 fcT =c feT , C H = (-lfic ki (6) 

where N-\ is the number operator of the up-spin fermions 
in the system. Apparently, these operators satisfy the 
conventional fermionic anticommutation relations 

i C ktT' = Skk', 

(7) 

{Ck«,C k ,„} = {Cl,Cij = 0. 

However, the operators C k -\ and C k >i now commute with 
each other. Consequently, Hamiltonian (|^) can be rewrit- 
ten in the following direct-product form, 

H M =^2^k (h k] <g> Ii + If <g> h k ^j 

k, a 

e fc |, \e k ,\<uj D 

Using the pseudo-fermion operators, we can also write 
the wavefunction of a state with iVj spin-up particles and 
Ni spin-down particles in the form, 

D T D l 

*=EE^™®xi. (9) 

m = l n— 1 
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In Eq. (||), Xs is a state vector denned by 

tf = 3L---<5L,|o) (io) 

where (fci, • • • , fc/,), £ = iVf, for er =f; L = Ni, for <r =j, 
denote the admissible states occupied by the fermions 
with spin er. Apparently, the entire set {xj} gives a nat- 
ural basis for V a (L), the subspace of L fermions with spin 
a, whose dimension is denoted by D a (L). 

For the case £)f = D±, it can be shown (see below) 
that the coefficient matrix W in (^) can be brought into 
a diagonal form through (separate) unitary transforma- 
tions in the V[(N-[) and Vj,(iVj_) subspaces, respectively. 
To generalize this property to the case ^ D± or equiv- 
alently iVf ^ N\ (which is the case for an odd number 
of electrons in the grain), we consider a larger subspace 
U X = V T (iV T )eF T (iV;) and^ x = ^(iV T )©n(^ilforthe 
spin-up and spin-down particles, respectively |17|1 . This 
way the number of basis D = D-f + for spin-up and 
spin-down states becomes the same. The original wave- 
function (^|) can now be written in the form, 

* = WmnXi ® XL (11) 

m. n 

where 

is a D x D square matrix. 

Let us now introduce the polar factorization lemma in 
matrix theory. 

Lemma: Let A be an arbitrary (not necessarily Hcr- 
mitian) n x n matrix. Then, there are two n x n unitary 
matrices U, V and an n x n diagonal semi-positive defi- 
nite matrix R such that 

A = URV, 

r mm i = r m 8 m m' and r m > 0, m = 1, • • • , n. (13) 

The proof of this lemma can be found in any standard 
textbook of matrix theory |l8| and is also presented in 
the appendix of Ref . |l7|] . 

Applying this lemma to W — URV, we can rewrite 
©as 

D D 

* = Yj WmnXi ® Xl 

jn=l n— 1 
D D 

= Y,Y,( URV ^mnXl®X l n 
■m— 1 7i—l 
D 

= $>*a/®^, (w) 
i=i 

where {77} is a set of nonnegative numbers and 



D D 
m— 1 n— 1 

Since U and 1/ are unitary, {ipj} and {w} are also 
orthonormal bases in subspaces TL-\ and TLi , respectively. 
Furthermore, since \P is an eigenvector of iV-p and iV^, the 
following conditions should hold 

D 

(* | 7V T | *) - X] I # I V>l> = ^T. (16) 

z=i 

and 

D 

{V\N L \*)=Y,ri(<t>i \N\4>,) = N U (17) 
where we have used the normalization condition, 

D 

(*i*) = E r ? = L ( 18 ) 

In both Eqs. @ and @, the spin indices are dropped 
in the sums, because, in each equation, only one species 
of spin is involved. 

It is useful to note that, although the transformations 
U and V may mix states with different number of parti- 
cles, those states ipj and tj)\ which actually appear in the 
sum ([m]), i.e., those correspond to ri > 0, are eigenstates 
of the number operators and iVj. , respectively. Thus, 

N | tPi) = N-\ | N\(t>l)=N L \ dpi), for n > 0. 

(19) 

The proof is left as an exercise for the reader. 



B. Parity parameter in the ground state 

With these preparations, we now summarize our main 
result in the following theorem. 

Theorem: Let E${M) be the ground state energy of 
Hamiltonian (|^) with M particles and Ap be the parity 
effect parameter defined in Eq. ([!]). Then, Ap is strictly 
positive for any integer N subject to N < Ne- 

Proof: To prove this theorem, let us now consider the 
ground state of a system with 2N + 1 particles and a 
total spin S. Since S+ and S- commute with Hm, we 
can specialize on the case S z — | (iVj — iVjJ = |, i.e., 
A T = N+l and N L = N. Let | * (2A+1)} be the ground 
state wavefunction with the above properties. Using the 
representation (Q for | ^> (2N + 1)) and the Hamilto- 
nian (||) , we can write the ground state energy of 2N + 1 
particles in the form, 
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£ (27V+1) 

(*o(2JV + l) | H M | tf (27V + l)) 
f>f [(V/|f T |^> + (^|n|^ 



z=i 



> 



+ 



- E r ixn a (< I Qt(M') I ^>(< I Qi(M') I <>, 

(20) 

where 

fa =J2 £ kn k a, Qa(k, W) = V 'g(K k^C^Ck'a (21) 
k 

By applying inequality | ab |< i(| a | 2 + | 6 | 2 ) to each 
term in the second sum ol Eq. (GO) and dropping the spin 
indices, we obtain 

E {2N+ 1) 

-]Tr 2 ^ |f|Vi>+(^ |T|Vi) 
1 D 

-J^rf [(4>i |T|0,) + fa, |f |^> 

- ^ E I Q(fc,fc') | V^>(^ 2 fOCMOUO 
-■ J n.n,^ | | <fe)<<fe | Q(k,k') | <^). 

(22) 

The right-hand-side of the above inequality is identi- 
fied with the sum of the expectation values of Hm in the 
two states, 

D D 

i=i i=i 

where tpi and </>/ are the complex conjugates of ipi and </>/, 
respectively. To see this, we note that T a is hermitian 
and {Q a (k, k')} are realm, the representation chosen. A 
straightforward substitution establishes that inequality 
( p2| ) is equivalent to 

E {2N + 1) > | H M | *0 + i(* 2 | tf M | tf 2 ). 

(24) 

^FromEqs. (||), we see that, by construction, and 
"J 2 are wavefunctions in the subspaces V{N^ — TVj — 
N + 1) and V(N^ = N± = TV), respectively. It is easy to 
verify that they are also normalized, 



(tfl | *l) = (* 2 I *2> =E r ' 2 = L 



(25) 



Therefore, we may regard and Vf^ as variational wave- 
functions for a system of 27V + 2 and 27V particles, re- 
spectively. By the variational principle, we obtain, 

£ (27V+1) > i^o(27V + 2) + is (27V). (26) 

In other words, Ap > 0. To finish our proof, we need to 
show that inequ ality ( |2^ ) is actually strict. 

Suppose Eq. (|22|) holds as an equality, we must then 
have, 



(i> h \clc k , \i> h ) = (<f>i 2 \qp w \<t> h ), 



(27) 



for all terms in (|2fj| ) with n t ^ 0, n 2 / and g(k, k') > 0. 
To show that this is not possible, let us consider a sub- 
set of terms with k = k' and l\ = Now, CtCk' = 
C[C k = n h . Applying Eqs. @ and @ to * (27V+1), 



we obtain 



D 



^^ r ^|n t |W=«ll 



(28) 



k 1=1 



and 



D 



(29) 



k 1=1 



1=1 



(Note that the expectation value of h k is real in any given 
state.) Since there are equal number of terms in the two 
sums, one must have at least one pair which violate Eq. 
(p7|). Consequently, when g(k,k') > for all pairs of 
admissible states with k = k', inequality ( prf ) must be 
strict, i.e., Ap is strictly larger than zero. 
Our proof is accomplished. QED. 



C. Lower bound 

In the above proof, we only showed that the parity 
parameter Ap is positive and nonzero. Actually, when 
(?(fc, k') — g > is a constant for pairs of the admissi- 
ble states, we can further derive a positive lower bound 
for Ap. For this purpose, we notice that inequality ( p2] ) 
and the ensuing inequality (|26j) are derived by replac- 
ing terms of the type ab in the second sum in Eq. ( |20| ) 
with (|a| 2 + |6| 2 )/2. Taking into account the fact that 
E (2N + 1) is real, we can express the error caused by 
this procedure in the form | a — b* | 2 /2 for each term 
in the sum. For simplicity, let us only estimate the error 
produced by handling terms with k = k' and l\ = I2 = I. 
For such terms, C\C k = n k . Combining Eqs. ©, (||) 
and (]24|), we obtain the following inequality, 

D 

k 1=1 
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Furthermore, by using the Cauchy-Schwarz inequality 

En a nbnf < (En I <*» I") ' (En I K P) , We obtain 



An > Z. 



g |E fc E/ r / 2 ((ipi I "fc I V'i) - {<t>i I I <&))| 



Efc Ez ' / 



iFrom Eqs. ©, © and (§|), we obtain finally, 



A P > 



2N f 



(31) 



(32) 



Here again Ap is the number of single-particle admissible 
states in the relevant energy range. 



IV. DISCUSSION AND CONCLUSIONS 

In this paper, by applying a generalized version of 
Lieb's spin-reflection positivity technique, we established 
the existence of a nonvanishing parity parameter Ap in 
an ultra-small superconducting grain and derived a pos- 
itive lower bound for Ap in a special case. Our study 
complements previous analytical work in the weak cou- 
pling regime 11 and numerical investigations [ p~T|JT2) . The 
bound ( p2| ) may seem to be too crude for practical pur- 
poses. However, we emphasize that it is rigorous and its 
validity does not require the coupling constant g to be 
small. 

An interesting open question is how pair-breaking scat- 
tering processes influence Ap in the ultra-small grain 
limit. This question has not been addressed explic- 
itly in previous theoretical studies |^-|l3]| based on the 
BCS Hamiltonian (|J), which does not contain any pair- 
breaking interaction. In contrast, experimental work 
by Ralph, Black, and Tinkham jl],|| suggests that pair- 
breaking due to spin-orbit interactions may be quite 
strong in real systems. Since our method relies only on 
certain global symmetries of the Hamiltonian, qualitative 
answers to the above question can be obtained in some 
cases. To be more precise, let us consider separately the 
following two situations. 

(I) Scattering processes that preserve spin symmetry. 
In such a process, the spin of the scattered electron is not 
flipped. A familiar example is the scattering of electrons 
by a nonmagnetic impurity described by the Hamilto- 
nian, 



H' = 



(t(k u k 2 )c{ ua c k2t „ + C.c)j , (33) 



where the scattering matrix (t(fci, fe)) is Hermitian. Ap- 
parently, H' is pair-breaking but it preserves the spin of 
the scattered electrons. We notice that the combined 
Hamiltonian H = Hm + H' can be formally thought to 
be an "extended negative-U Hubbard model" defined on 
a lattice in the single- particle eigenstate space. (In the 
literature of the strongly-correlated fermion models, such 



a lattice is called a "superlattice" [ |19[]20| ] ) . Such a model 
can be treated in a similar way as in Ref . fl7|| , where the 
binding energy of the negative-?/ Hubbard model on a 
real-space lattice was discussed. By following the proof 
of theorem I in Ref. [|l7j and the proof of the above in- 
equality (|32]), without further ado, one easily finds that 
the main result of this paper is still valid. Thus, the 
scattering of electrons by nonmagnetic impurities does 
not destroy superconductivity in a small superconduct- 
ing grain. 

(II) Scattering processes that do not preserve spin sym- 
metry. In such a process, the Cooper pair is broken and 
the spin of the scattered electron is also flipped by the 
interaction. Examples include scattering by a magnetic 
impurity and spin-orbit interactions. While the former 
process definitely suppresses superconductivity in the su- 
perconducting grain, the origin and effect of the latter 
interaction have not been fully understood [QJ^]. Due 
to the spin-flipping process, the total electron spin of 
the system is no longer a good quantum number. Con- 
sequently, our approach based on Lieb's spin-reflection 
positivity is not applicable to these systems. To analyze 
the effects of such interactions on the superconducting 
grain in a mathematically rigorous way, more sophisti- 
cated techniques are needed. Further research on this 
issue is currently under way. 
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